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^ ! Abstract 
(N 

' Single scale quantities, as anomalous dimensions and hard scattering cross sections, in 

■ renormalizable Quantum Field Theories are found to obey difference equations of finite 

order in Mellin space. It is often easier to calculate fixed moments for these quantities 
compared to a direct attempt to derive them in terms of harmonic sums and their gen- 
^ ' eralizations involving the Mellin parameter N. Starting from a sufficiently large number 

of given moments, we establish linear recurrence relations of lowest possible order with 
polynomial coefficients of usually high degree. Then these recurrence equations are solved 
in terms of d'Alembertian solutions where the involved nested sums are represented in op- 
timal nested depth. Given this representation, it is then an easy task to express the result 
in terms of harmonic sums. In this process we compactify the result such that no algebraic 
relations occur among the sums involved. We demonstrate the method for the QCD un- 
polarized anomalous dimensions and massless Wilson coefficients to 3-loop order treating 
the contributions for individual color coefficients. For the most complicated subproblem 
5114 moments were needed in order to produce a recurrence of order 35 whose coefficients 
have degrees up to 938. About four months of CPU time were needed to establish and 
solve the recurrences for the anomalous dimensions and Wilson coefficients on a 2 GHz 
machine requiring less than 10 GB of memory. No algorithm is known yet to provide such 
a high number of moments for 3-loop quantities. Yet the method presented shows that it 
is possible to establish and solve recurrences of rather large order and degree, occurring in 
physics problems, uniquely, fast and reliably with computer algebra. 
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1 Introduction 



Precision predictions for observables in Elementary Particle Physics require the calculation of 
the corresponding Feynman diagrams, the number of which grows fast with the order in the 
coupling constant being considered. According to the relevant number of different ratios of 
Lorentz invariants or scales involved one may group these observables into 0-scale, 1-scale, 2- 
scale etc processes. In renormahzable Quantum Field Theories the radiative corrections to 
the couplings, masses and external fields are examples for 0-scale quantities [1]. Anomalous 
dimensions and hard scattering cross sections, as the Wilson coefficients for light and heavy 
flavors (for ^ m|^) in deeply inelastic scattering, are single scale quantities, cf. [2-6]. Also 
the sub-system cross sections for the Drell-Yan process and the cross section for hadronic Higgs- 
boson production in the heavy mass limit for the top-quark belong to this class. Mellin moments 
for single scale quantities /(x), 

M[/(x)](Ar)= dxx"" j{x) (1) 
Jo 

are 0-scale quantities again for G N [7-10]. Here x usually denotes a fraction of Lorcntz- 
invariants the support of which is or can be extended to [0, 1]. In the lower order in perturbation 
theory 0-scale quantities can be expressed as linear combinations of specific numbers over Q 
which are multiple C- values [11], 

Cai....,an = E «^ e Z \ {0} 

fc=l 

at the beginning, with possible extensions in higher orders, which occur in both massive and 
massless calculations [12]. The 1-scale quantities can be expressed in terms of finite harmonic 
sums [13, 14] 

^ ■ I \k 
fc=l 

and rational functions of the Mellin variable A^ at lower orders in perturbation theory. At higher 
orders one expects to find generalizations of harmonic sums. Much less is known on the function- 
spaces spanning 2- and higher scale processes. The Mellin-transformation (1) is empirically found 
to yield considerable structural simplifications of 1-scale processes, cf. [15]. In massless processes 
this is partly due to the factorization properties, but it seems to hold to an even wider extent. 
Corresponding diagonalizations for processes with a higher number of scales depend on their 
respective main symmetries, which may not even be fixed by just the number of scales. 

In the present paper we study single scale processes and represent them in Mellin space. In 
order to apply our method under consideration, we shall assume the case that M[/(a;)](A^) can 
be found as the solution of a linear recurrence equation 

ao(iV)F(iV) + a^{N)F{N + 1) + • • • + ai{N)F{N + /) = 0, (2) 

with polynomial coefficients ak{N). 

There is no general proof that the /c-loop contributions to a 1-scalc observable have to obey 
such a recurrence. On the other hand, it is known that all single scale processes having been 
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calculated so far do, cf. [2-5, 15]. This is due to the fact that the corresponding observables are 
found as linear combinations of nested harmonic sums. The single harmonic sums obey 



Exploiting holonomic closure properties [16] one obtains higher order difference equations for 
polynomial expressions in terms of nested harmonic sums. 

If a suitably large number of moments NL[f[x)]{N) is known, then a recurrence of the form (2) 
can be found automatically, see Section 2. Once a recurrence of some order I is found, this 
recurrence together with the first I moments specifies uniquely all the moments lSA[f{x)][N) 
for nonnegative integers N. Finally, we activate the summation package Sigma [17] and solve 
the recurrence (2) in terms of generalized harmonic sums. In particular, using the underlying 
summation theory of IIE-difference fields [18-20] or exploiting the algebraic relations [21], a 
closed form for 'M[f(x)]{N) in terms of an algebraically independent basis of harmonic sums 
can be computed. 

Wc emphasize that Eq. (2) covers a much wider class in which more general recurrent quanti- 
ties can represent the corresponding observables. In particular, our general recurrence solver for 
d'Alembertian solutions [22-24] finds any solution that can be expressed in terms of indefinite 
nested sums and products. In even higher order or massive calculations further functions may 
contribute, which could be only found in this way. 

The Mellin-moments of the unpolarized 3-loop splitting functions and Wilson coefficients 
for deep-inelastic scattering are more easily calculated [7-10] than the complete expressions, 
cf. [3-5]. In the present paper we investigate whether the exact formulae up to the unpolarized 
3-loop anomalous dimensions and Wilson coefficients [3-5] can be found establishing and solving 
difference equations (2) for the Mellin moments of these quantities, without further assumptions.^ 

We consider the various color contributions to these quantities separately and try to find the 
complete result from a minimal number of moments. As input we apply the moments calculated 
from the exact solution [3-5] . 

The paper is organized as follows. In Section 2 we describe how the difference equations of 
the form (2) are found by just using a finite number of starting points of F{N). In Section 3 the 
algorithms are outlined that can solve these recurrences in the setting of difference fields. They 
lead directly to the corresponding mathematical structures. These are nested harmonic sums 
in the present case. In course of the solution we compactify the results applying the algebraic 
relations to the harmonic sums [21]. ^ The results are discussed in Section 4. Our method applies 
in the same way to all other single scale processes of similar complexity, cf. [6, 15]. Section 5 
contains the conclusions. In the appendix we present a compactified form of the non-singlet 
3-loop anomalous dimensions, which is automatically provided in the formalism by Sigma. The 
corresponding expressions for the other anomalous dimensions and Wilson coefficients to 3-loop 
order are presented in Mathematica and FORM codes attached. 

"'^Approximate reconstruction methods based on special ansatzes were discussed in the literature e.g. in [7,25] 
to obtain first numerical estimates from a low number of moments, see also [26]. We also remind that the 
description of QCD-cvolution relating fixed integer moment information to orthogonal polynomials is an old 
topic [27]; see also [28]. 

^Further compactifications can be obtained using the structural relations, cf. [29,30]. 
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2 Finding a Recurrence Equation 



Suppose we are given a finite array of rational numbers, 

qi, ?2, Qk , 

which are the first terms of a certain infinite sequence F{N), i.e., = qi, F[2) = ^2, etc. Let 
us assume that F{N) satisfies a recurrence of type 

I d 

Y.{^c,,,N^)F{N + k)^Q, (3) 

A;=0 1=0 

which wc would like to deduce from the given numbers qi (i = 1, . . . , K). In a strict sense, this 
is not possible without knowing how the sequence continues for N > K. One thing we can do 
is to determine the recurrence equations satisfied by the data we are given. Any recurrence for 
F{N) must certainly be among those. 

To find the recurrence equations of F{N) valid for the first terms, the simplest way to 
proceed is by making an ansatz with undetermined coefficients. Let us fix an order / G N and 
a degree G N and consider the generic recurrence (3), where the q^/j are indcterminates. For 
each specific choice N = 1,2, . . . , K — I, we can evaluate the ansatz, because we know all the 
values of F{N + k) in this range, and we obtain a system oi K — I homogeneous linear equations 
for {I + l){d + 1) unknowns Cjj. 

li K — I > {I + l){d + 1), this system is under- determined and is thus guaranteed to have 
nontrivial solutions. All these solutions will be valid recurrences for F{N) for N = 1, . . . , K — I, 
but they will most typically fail to hold beyond. If, on the other hand, K — I < (/ + l)((i + 1), 
then the system is overdetermined and nontrivial solutions are not to be expected. But at 
least recurrence equations valid for all N , if there are any, must appear among the solutions. 
We therefore expect in this case that the solution set will precisely consist of the recurrences 
of F{N) of order / and degree d valid for all A^. 

As an example, let us consider the contribution to the gluon splitting function oc Ca at 
leading order, Pgg\N). The first 20 terms, starting with A^ 3, of the sequence F{N) are 

14 21 181 83 4129 319 26186 18421 752327 71203 811637 128911 29321129 

5 ' 5 ' 35 ' 14' 630 ' 45 ' 3465 ' 2310 ' 90090 ' 8190 ' 90090 ' 13860 ' 3063060 ' 

2508266 292886261 7045513 611259269 1561447 4862237357 988808455 

255255 ' 29099070 ' 684684 ' 58198140 ' 145860 ' 446185740 ' 89237148 ' 

Making an ansatz for a recurrence of order 3 with polynomial coefficients of degree 3 leads to 
an overdetermined homogeneous linear system with 16 unknowns and 17 equations. Despite 
of being overdetermined and dense, this system has two linearly independent solutions. Using 
bounds for the absolute value of determinants depending on the size of a matrix and the bit 
size of its coefficients, one can very roughly estimate the probability for this to happen "by 
coincidence" to about 10"^^. And in fact, it did not happen by coincidence. The solutions to 
the system correspond to the two recurrence equations 

{7N^ + nSN^ + 494A^ + 592)F{N) - {12N^ + 233N^ + 1289Ar + 2156)F(A^ + 1) 

+ {3N^ + IISN^ + 1021A^ + 2A76)F{N + 2) + {2N^ + 2N^ - 226N - 912)F{N + 3) = (4) 

and 

{AN^ + QAN^ + 278N + 332)F{N) - {7N^ + 13AN^ + 735N + 1222)F{N + 1) 

+ {2N^ + 71N^ + 595iV + 1418)F(A^ + 2) + (A^^ - A^^ - 138iV - 528)F(A^ + 3) = 0, (5) 
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which both are vahd for all > 1. If we had found that the hnear system did not have a 
nontrivial solution, then we could have concluded that the sequence F{N) would definitely (i.e. 
without any uncertainty) not satisfy a recurrence of order 3 and degree 3. It might then still 
have satisfied recurrences with larger order or degree, but more terms of the sequence had to be 
known for detecting those. 

The method of determining (potential) recurrence equations for sequences as just described 
is not new. It is known to the experimental mathematics community as automated guessing 
and is frequently applied in the study of combinatorial sequences. Standard software packages 
for generating functions such as gfun [16] for Maple or GeneratingFunctions.m [31] for Mathe- 
matica provide functions which take as input a finite array of numbers, thought of as the first 
terms of some infinite sequence, and produce as output recurrence equations that are, with high 
probabihty, satisfied by the infinite sequence. 

These packages apply the method described above more or less literally, and this is perfectly 
sufficient for small examples. But if thousands of terms of a sequence are needed, there is no 
way to get the linear systems solved using rational number arithmetic. Even worse, already for 
medium sized problems from our collection, the size of the linear system exceeds by far typical 
memory capacities of 16-64Gb. For the big problem C^^^ „3 (A^), it would require approximately 
11Tb of memory to represent the corresponding linear system explicitly. It is thus evident 
that computations with rational numbers are not feasible. Instead, we use arithmetic in finite 
fields together with Chinese remaindering and rational reconstruction [32-34]. Modulo a word 
size prime, the size of the biggest systems reduces to a few Gb, a size which easily fits on our 
architecture. And modulo a word size prime, such a system can be solved within no more than 
a few hours of computation time by Mathematica. 

The modular results for several distinct primes pi,P2, ■ ■ ■ can be combined by Chinese re- 
maindering to a modular result whose coefficients are correct modulo the product pip2 .... If 
the bit size of this product exceeds twice the maximum bit size appearing in the rational solu- 
tion, then the exact rational number coefficients can be recovered from the modular images by 
rational reconstruction [32-34]. The number of primes needed (and thus the overall runtime) is 
therefore proportional to the bit size of the coefficients in the final output. 

The final output is a recurrence equation for F{N). But the recurrence equation satisfied 
by a sequence F{N) is not unique: if a sequence satisfies a recurrence equation at all, then 
it satisfies a variety of linearly independent recurrence equations. The bit size of the rational 
number coefficients in these recurrence equations may vary dramatically. In order to minimize 
the number of primes needed for the computation of the rational numbers in the recurrence, it 
seems preferable to compute on a recurrence whose coefficients are as small as possible in terms 
of bit size. According to our experience, this recurrence happens to be the (unique) recurrence 
whose order I is minimal among all the recurrence equations satisfied by F{N). We have no 
explanation for this, but it seems to be a general phenomenon, as it can also be observed in 
certain combinatorial applications [35]. 

Also the number of unknowns for the linear system may vary dramatically among the possible 
recurrence equations for F{N), and it seems preferable to compute on a recurrence where the 
number of unknowns is as small as possible. Small linear systems are not only preferable because 
of efficiency, but also because the number of unknowns in the linear system determines the 
number of initial terms qi that have to be known a priori in order to detect the recurrence. 
According to our experience, the size of the linear system is minimized when the order / and the 
degree d are approximately balanced. 

Unfortunately, it seems that the recurrence with minimal (in terms of bit size) rational 
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number coefficients lias tfie maximal number of unknowns in tfie corresponding linear system, and 
vice versa. But there is a way to combine the advantages of both at a reasonable computational 
cost. Consider the two recurrence equations (4) and (5) from the example of the gluon-gluon 
splitting function at leading order, Pgg Q{N), quoted above. A recurrence of smaller order can 
be obtained from these by multiplying (4) by (A^^ - 9A^ - 66) and (5) by {2N'^ - UN - 114), 
and then subtracting the results. The choice of the multipliers is such that the coefficient of 
F{N + 3) in the difference cancels: we obtain 

(iV^ + 22N^ + 189iV^ + 788iV^ + 1592iV + 1224)F(iV) 

- (2iV^ + 45iV^ + 396iV'^ + ITOliV^ + 3580A^ + 2988)F(iV + 1) 

+ {N^ + 23N'^ + 207N^ + 913iV^ + 1988Ar + 1764)F(Ar + 2) = 0. 

The calculation just performed can be recognized as the first step in a difference operator ver- 
sion of the Euclidean algorithm [36]. Applied to two recurrence equations satisfied by a se- 
quence F[N), this algorithm yields their "greatest common (right) divisor", which is, with high 
probability, the minimal order recurrence satisfied by F{N). In our example, the algorithm 
terminates in the next step, and indeed the sequence F{N) of Pggfi{N) does not satisfy a recur- 
rence of order less than two. Note that the linear system for finding the second order recurrence 
directly would have involved (5-|- 1)(2-|- 1) = 18 unknowns instead of the 16 unknowns we needed 
for finding the third order recurrences. For the big problem C^^^^-i{N), a direct computation 
would require 33804 unknowns instead of the 5022 wc actually used. We combine the advantage 
of a small linear system with the advantage of small coefficients in the output as follows. We 
first compute for several word size primes the solutions of a small linear system, but then instead 
of applying rational reconstruction to those, wc compute, for each prime independently, their 
greatest common right divisor modulo this prime. We then apply rational reconstruction to 
recover the rational number coefficients of those. 

In summary, we used the following procedure for finding the recurrence equations. 

1. Choose a word size prime p. 

2. Choose some bounds I and d and make an ansatz for a recurrence of order I and degree d. 
The linear system is constructed and solved modulo p only. 

3. If there are no solutions, repeat step 2 with increased bounds I and d. 

4. If there are solutions modulo p, compute their greatest common right divisor modulo p by 
the Euclidean algorithm for difference operators. 

5. Repeat steps 1-4 until Chinese remaindering and rational reconstruction applied to the 
greatest common right divisors for the various primes yields a recurrence that matches the 
given data qi, q2, ■ ■ ■ , Qk- 

6. Return the reconstructed recurrence as the final result. 

For the big problem C2^^^^3(A^), most of the computation time (about 53%) was spent in 

step 4. Solving the modular linear systems consumed about 28% of the time, and Chinese 
remaindering and rational reconstruction took about 18% of the time. The memory bottleneck 
is in step 2 where the linear system is constructed. The memory requirements for the other 
steps, if implemented well, are negligible. 
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For problems that are even bigger than those we considered, further improvements to the 
procedure are conceivable. First, there are asymptotically fast special purpose algorithms for 
step 2 available [37,38]. These algorithms outperform the naive linear system approach we are 
taking for problem sizes where fast polynomial multiplication algorithms outperform classical 
algorithms. It is likely that their use would have already been beneficial for some of our problems. 
Second, a gain in efficiency might result from running the procedure on a different platform. We 
have done all our computations within Mathematica 6, but we expect that in particular step 4 
might considerably benefit from a rcimplementation in a computer algebra system providing 
high-performance polynomial arithmetic. Mathematica 's modular arithmetic, on the other hand, 
appears to be quite competitive. Third, it might be worthwhile to run parts of the procedure in 
parallel. In particular, computations for distinct primes are completely independent from each 
other and can be done on different processors without any communication overhead. Observe 
that these steps dominate the runtime. 



After having obtained difference equations of high order and degree we will now discuss general, 
efficient algorithms by which these equations can be solved. Given a recurrence relation 



of order I, find all its solutions that can be expressed in terms of indefinite nested sums and 
products. Such solutions are also called d'Alembertian solutions [22-24], they form a subclass 
of Liouvillian solutions [39]. Note that such solutions cover as special cases, e.g., harmonic 
sums [13,14] or generalized nested harmonic sums [40]. 
The solution to this problem consists of two parts. 

1. First, compute all d'Alembertian solutions by factoring the recurrence as much as possible 
into linear right factors. Then each linear factor contributes to one extra solution. To be 
more precise, the ith factor yields a nested sum expression of depth i — 1. 

2. Second, simplify these nested sum solutions to closed form expressions, e.g., in terms of 
harmonic sums, that can be processed further in practical problem solving. 

In general, the package Sigma [17] can solve these problems in the setting of IlS-difference 
fields [18,41]. This means that the coefficients ao{N), . . . ,ai{N) and the inhomogeneous part 
q{N) of (6) can be given as polynomial expressions in terms of indefinite nested sums and 
products. 

For simplicity, we restrict ourself to the situation that the given coefficients ao{N), . . . , ai{N) 
are polynomials in N and that the inhomogeneous part q{N) is zero. In other words, we assume 
that we are given a recurrence of the form (2) or (3) that is produced, e.g., by the method 
described in the previous section. 

3.1 Finding all d'Alembertian solutions 

Subsequently, we present algorithms that find all d'Alembertian solutions of (2). Equivalently, 
we can say that we look for all d'Alembertian sequences which are annihilated by the linear 



3 Solving the Recurrence Equations 



ao{N)F{N) + ai{N)F{N + 1) + • • • + ai{N)F{N + I) ^ q{N) 



(6) 



operator 



L := ao{N) + ai{N)S + ■■■ + ai{N)S\ 



(7) 
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which is understood to act on a sequence F{N) via 

(L • F){N) := ao{N)F{N) + ai{N)F{N +!) + ••• + ai{N)F{N + I). 
We start as follows. 

Step 1: Finding a product solution. First, we look for a solution of (7) which is of the form 

N 



To{N)=\{r{i) (8) 



for some rational function r{i) in i. In Sigma this task can be carried out by executing a 
generalized version of algorithm [42] that works in general IIE-difference fields; for an alternative 
algorithm to find such hypergeometric terms we refer to [43] . 

If there does not exist such a product solution (8). then there is no d'Alembertian solution 
at all; see, e.g. [24, Theorem 4.5.5]. In this case, we just stop. Otherwise, we look for additional 
solutions as follows. 

Step 2: Splitting off a linear right factor. By dividing the operator (7) from the right with the 
operator 

we arrive at an operator 

L' := bo{N) + 6i(iV)S + • • • + 6i_i(iV)S'-i (10) 

of order I — 1 such that 

L = L'{S-r{N + 1)) 

= -r{N + l)bo{N) + {bo{N) - r{N + 2)b,{N))S + ■■■+ {bi.,{N) - r{N + l)bi{N))S\ 
i.e., S — r{N + 1) is a linear right factor of L. 

Step 3: Recursion. Now we continue by recursion and look for all d'Alembertian solutions for 

the operator L' with order / — 1. Note that after at most / — 1 steps we end up at a recurrence 
of order 1 whose d'Alembertian solution can be read off immediately. 

Step 4: Combining the solutions. If we do not find any d'Alembertian solution for L', we just 
return the solution (8) for L. 
Otherwise, let 

ti(iV),...,4(iV) (11) 

with 1 < k < I he the solutions of L' that we obtained after the recursion step. To this end, for 
1 < i < ^ define 

N . /• _ -|X 

for some properly chosen A > (i.e., To(i) is nonzero for all i with i > X). Then the final output 
of our algorithm is 

n{N),T,{N),...,n{N). (13) 
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The following remarks are in place. By construction all the elements from (13) are solutions 
of (7): for each 1 < j < k, 

(S-r(iV + 1)) . T,(iV) = To{N + 1) ^^f-ii - r(N + l)n(N) ^ 



Toil) 

^ riN + DUN) f ± ^ + ^f^^ ]-riN+ l)To(iV) f ^ = t,(iV) 




N 



i=\ 



and hence 



L-T,(iV) = L'-t,(iV)=0. 

But even more holds. The derived solutions (13) are linearly independent. In particular, any 
solution of L in terms of indefinite nested sums and products can be expressed as a linear 
combination of (13); see [39, Theorem 5.1] or [24, Proposition 4.5.2]. 

Summarizing, with the algorithm sketched above we can produce all d'Alembertian solutions 
of L, i.e., all solutions that are expressible in terms of indefinite nested sums and products. 

Wc emphasize that the expensive part of the sketched method is the computation of the 
product sohitions (8). The following improvements were crucial in order to solve the recurrences 
under consideration. 

Improvement 1. If one finds several product solutions, say Pi{N), . . . , Pu{N), one can produce 
immediately a recurrence L' like in (10), but with order l—u instead of order l—l. Moreover, given 
all d'Alembertian solutions of this operator L', one gets all the solutions of the recurrence (7) 
without any further computations; see [24, Theorem 4.5.6]. 

Improvement 2. For the problems under consideration, it turns out that it suffices to search for 
product solutions (8) that can be written in the form 

^ o. (14, 

for polynomials p{N) and q{N). Therefore, we used optimized solvers [44] of Sigma which 
generalize the algorithm presented in [45]. In addition, arithmetic in finite fields is exploited in 
order to determine the solutions (14) effectively. 

Improvement 3. In our applications, rather big factors from tj{i — 1) and T'o(i) cancel in the 
summand of (12); in particular, the usually irreducible factor p(i) from (14) (A^ substituted 

with i) cancels. Hence it pays off to compute directly the summand expression '■ Namely, 

instead of the operator (10) we continue with the operator 

ho{N) + r{N)h{N)S + ---+[^r{N + i) j hi^^{N)S'-\ (15) 

and look for all its d'Alembertian solutions t[{N), . . . , t'i^{N). Then by construction, the solutions 
of (6) can be given directly in the form 

N 

Tj{N)^To{N)J2m- (16) 

i=X 
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Example 1 As illustrative example we solve the difference equation for the CpN^-term of the 
unpolarized 3-loop splitting function Pgq^2{N), 

F{N) = Pgq^,^^.cAN) . 
Using the methods from the previous section, we generate the recurrence relation 

ao{N)F{N) + ai{N)F{N + 1) + a2{N)F{N + 2) + a3{N)F{N + 3) = , (17) 

with 

ao{N) = (1 - N)N{N + 1){N^ + 15N^ + mN^ + A85N^ + 1358N^ + 2216A^ + 1616), 
a^(N) = N{N + 1){3N^ + ASN^ + 366N^ + 17A0N^ + 5527N^ + 11576A^^ + 14652Ar + 8592), 
02 (AT) = -(AT + 1){3N^ + MN^ + A57N^ + 2441 A^^ + 9064A^^ + 23613A^^ 

+ 41180A^2 + 43172A^ + 20768), 
a3{N) ^{N + Af{N^ + 9N^ + 49A^^ + 179A^^ + 422A^2 + 588A^ + 368) . 

Given this recurrence, we produce its d'Alembertian solutions as follows. First, Sigma computes 
a rational solution, namely 

T (N) = ^" + ^ + ^ 
°^ ^ {N -l)N{N + iy 

Now we can divide (17) from the right by the operator S—To{N +1) /To{N) — S 



{N+2){N'^+N+2) ■ 

Then the resulting recurrence of the operator (15) is 

bo{N)G{N) + bi{N)G{N + 1) + b2{N)G{N + 2) = , (18) 

with 

bo{N) = {N + 1){N + 2){N^ -N + 2){N'^ + 9N^ + A9N^ + 179N^ + A22N^ + 5SSN + 368), 
b^{N) = -{N + 2){2N^ + 21 A^^ + 124A^^ + b3m^ + 1690A^^ + 3989 A^^ + 6712Ar3 ^ 7524A^^ 
+ 5232Ar + 2080), 

b2{N) = (A^2 + 5A^ + 8)(A^^ + 3A^^ + 19A^^ + 53A^^ + IMN'^ + 124A^ + 64) (AT + 3)^ . 

Next, we proceed recursively and can compute the rational solution 

+ AN^ + 13Ar2 + 22Ar + 16 



P\N) 



{N - 1)N{N + 1){N + 2) {N^ + 3Ar + 4) 



of (18). Thus we divide (18) by the factor S — ^p!^^^^ which leads to the first order recurrence 

co{N)H{N) + ci{N)H{N +1) = (19) 



with 



Co{N) = -{N + 1){N'^ + N + 2){N^ - 4N^ + 13N^ - UN + 8) 

X {N^ + 3N^ + 19Ar^ + 53N^ + lOAN^ + 12AN + 64) , 

ci{N) = {N + 2){N^-N + 2){N^ + AN^ + 13N^ + 22N + 16) 
X {N^ - 3N^ + 19iV^ - 13Ar3 + 44iV2 + 8iV + 8). 
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Here we can read off directly the solution 

{N^ -N + 2) {N^ - 3N^ + IQiV^ - 13N^ + AAN^ + 8N + 8) 



P"{N) = 



{N + 1) (iV4 + + AN + A) {N^ - 47V3 + i^p^2 _ ^4;^ + g) • 
Going back, we obtain besides ti{N) — P'{N) the solution 

N ^ / X / X 

f (AT) _ p'(]V\\^P"(l\ N^+7N^+4N+A (i^-j+2)(/-3/+19,-^-13/+44j^+8j+8) 

L2{1^ ) - -r {1\ ) 1- {J) - (jv+l)(Ar^-iV+2)(iV^+iV+2) (j+i)(j4+7j2+4j+4)(j-4_4j-a+i3j2_i4j+8) 

0/ (18). Hence by (16) we obtain, besides Tq{N) , the solutions 

^ i^ + 7i^ + Ai + A 



^ (i + l)(i2-i + 2)(i2 + i + 2) 



{N - 1)N{N + 1) 

C,-4 , 7,-2 , 4,- , 4^ V (i^-i+2)(i^-3i^+19i^-13j^+44i^+8j+8) (20) 
AT ^ 0+l)0-*+7j2+4j+4)0-4_4j3+i3j2_i4j+8) 



7^2(71) 



(Ar2 + AT + 2) y 1± 

^ ■fcl' (^ + l)(^2-2 + 2)(z2+Z + 2) 

■ {N -1)N{N + 1) 



for (17). S'mce all three solutions Tq{N), Ti{N) and T2{N) are linearly independent over, say, 
the complex numbers, any solution F : N —>■ C of (17) can be described as a linear combination 

F{N) = ciTo(7V) + c2Ti(7V) + csT^iN) 

for ci, C2, C3 e C. The initial values F{3) = F(4) = F(5) = ||| zm;j/y 

n,,2,iv|c.(A^) = -fTo{N) + f Ti(iV) - |r2(Ar). (21) 

For our concrete problems all the recurrences could be factored completely. Equivalently, for 
a recurrence of order d we found d linearly independent solutions Ti(N), . . . ,T(i(N) where the 
solution Tk with 1 < k < d can be given in the form 

T,(m- V ^'^''^ V ^^^'^^ V ^"^'''^ (22) 

<Vo(A^);^^ '^U^i)^^! Q2{i2) Qkiik) 
where for 1 < i < /c the Pj and are polynomials and Sj e {—1, 1}. 

Example 2 For ^/le CAC^Np-term of the 3-loop non-singlet splitting function P^S2 found 

a recurrence of order 7 which fills around five pages. The 7 linearly independent solutions can 
be computed within 10 seconds; the largest solution fills around three pages and has the form 

2^ nj^\ ^ n^( i\ ^ njh\ 2^ n.(i\ 2^ njy) 2^ nj <,\ ^ ' 



i=6 



where the irreducible polynomials Pi, P2, . . . , have the respective degrees 4, 8, 16, 28, 63, 69, 
and the denominators are of the form 

Qi{s) ={s-lfs{s + lf, 
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Q2{r) = (r^ - lOr^ + 29r^ - 34r + 12) (r^ - 6r^ + Sr^ - 2r - 2) , 

Q3(/) =(^8 - 24f + 215Z^ - 1017/^ + 2866Z'' - 4975/^ + 5146/^ - 2812Z + 576) 

X {l^ - 16f + 75f - 175/^ + 236/^ - 165/^ - 4/^ + 64/ - 24) , 
Q4{k) ={k - 4){k - 3f{k - 2f{k - lfk^{k^^ - 38k^ + 566k^ - 4628/c^ + 23621A;^ - 79466A;^ 

+ 178404A;^ - 261580A;^ + 235712A;^ - 114624A; + 21600) {k^° - 28k^ + 269A;^ - 1348A;^ 

+ 4091A;6 - 7768A;^ + 8451A;^ - 3560A;^ - 1612A;2 + 1872A; - 432) , 
Q^{j) =2/0 _ 795/9 ^ 40760i^^ - 1036641 j^'^ + 16752826^^^ - 191239786/^ + 1632641752/'' 

- 10786299042^^3 + 56334695030/^ - 235648109263^^^ + 795075807544j^° 

- 2168602473357^'^ +4771126881598/ - 8409573468828/ + 11731291260824^^ 

- 12705852943232 + 10375981856560/ - 6104512549760^3 + 2399836168064^2 

- 547585520256j + 51445094400, 

Qe(i) ={i - - + I)(l6z33 - 7192z32 + 673840?^! - 33108234?^° + 1069628658^^^^ 

+ • • • + 162245083333715039232Z - 11706508031797555200) (16^33 _ 6664z32 + 452144^3^ 

- 15699130^30 + ■ ■ ■ + 6071537402380800?^ - 670382971978752? + 32623028121600) . 

Example 3 The solution of the recurrence for the Cp-contribution to the unpolarized 3-loop 
Wilson coefficient for deeply inelastic scattering, c!.^^ ^^{N), constituted the hardest problem to 
solve. We obtained a recurrence of order 35. Then our solver ran 25 hours and used 3 GB of 
memory to derive the 35 linearly independent solutions. In total, we needed only J^IS instead 
of Yl^to ^ = 595 summation quantifiers in order to represent those solutions. This is possible 
due to the Improvement 1. For each of the summands around 20 MB of memory were used. In 
particular, in the summands the denominators have irreducible factors up to degree 1000; the 
integer coefficients of the polynomials were up to 700 decimal digits long. 

3.2 Simplification of d'Alembertian solutions 

We consider the following problem: Given indefinite nested sum and product expressions, e.g., 
expressions of the form (22), find an alternative sum representation with the following properties: 

1. All the involved sums are algebraically independent with each other. 

2. The nested depth of the sum expressions is minimal. 

3. In the summands the degree of the denominators is minimal. 

4. The sums should be tuned in such a way that algorithms can perform this simplification 
as efficiently as possible. 

In principal, this problem can be solved with Karr's summation algorithm [18] based on IIE- 
difference fields, if one knows explicitly the sum elements in which, e.g., the expression (22) should 
be expressed. For small examples such optimal sums with properties 1-3 from above might be 
guessed. In particular, if one has additional knowledge about the objects under consideration, 
a good sum representation might be known a priory. But if such additional knowledge is not 
available, Karr's algorithm is not applicable. 

In order to overcome this restriction, the fourth named author has refined Karr's IIE-theory 
for symbolic summation [20,46]. As a consequence, we can determine completely automatically 
such sum representations with the properties 1-4 from above; see [19,47]. 
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Example 4 With Sigma we find the depth- optimal representation 

4:{N^ + N + 2) 8 (8A^=^ + 13A^2 + 27A^ + 16) ^ 1 

~3(A^ - 1)A^(A^ + 1) [j^ 1 J ^ 9(A^ - 1)A^(A^ + 1)2 

8 {AN^ + AN^ + 23A^2 + 2bN + 8) 4 (A^^ ^ + 2) ^1 
9(A^ - 1)A^(A^ + 1)3 3(A^ - 1)A^(A^ + 1) ^ 

o/ (23) where the sums are given in (20) . We can read off the harmonic sum representation 
4(Ar2 + Ar + 2) ^^^^^^^ ^ 8(87V3^137V2 + 277V + 16) ^^^^^^^ 



3(A^ - l)A^(Ar + 1) ' ' ^{N -l)N{N + lf 

8 (4A^4 + AN^ + 23A^2 + 2bN + 8) 4 (A^^ ^ + 2) 



9(A^ - 1)A^(A^ + 1)3 3(A^ - l)Ar(A^ + 1) 



S2{N). 



Example 5 The sum expression for Pns,2,CaCfNf containing in particular the 7-nested sum 
(21) can be simplified with Sigma to the depth-optimal representation 

_ 2 (1086A^^ + 3258Ar6 + 2129A^^ - 288A^^ - 67N^ - 206N^ - 156N + 144) 
PnS,2,CaCfNf — 27N'^{N + 1)3 

32 (8Ar4 + 33A^3 ^ 53^2 ^ 25Ar + 3) , 16 1 32 A (-1)^ 



9Ar(Ar + i)4 ^ ^ ' 3 ■ 3 ^ ^4 

16 (lOAT^ + lOAT + 3) A (-1)' 1336 ^1 64 (8Ar2 + 8Ar + 3) A (-1)^ 
9iV(7VTT) ^ i3 ^ j-^ P 97V(7VTT) ^ i 

16 {AN^ + 88N^ + 314A^4 + 412A^3 ^ 201Ar2 + 16A^ - 12) A (-1)^ 
^ 9A^2(7v + l)2(Ar + 2)2 ^ ^2 

/ 8(14Ar2 + 14Ar + 3) 16f l\f 1 64 f F , S=i^ 
^1 3Ar(Ar + l) 3 z ) i3 ^ 3 i ^ £f (^ + 2)^ 



^=l / i=l 2=1 2=1 



32 {22N^ + 22N - 3) A E]=i V f 32 {2N^ + AN + 1) 

9A^(Ar + l) ^ i + 2 ^ i y \^ 3(iV+l)3 ^ ^ 

4 (65Ar6 ^ ;^95^5 ^ ;^95^4 ^ ;^37^3 ^ 35^2 ^ 35^ ^ ;^8) 33 ^ 128 A (-1)^ 



27Ar3(7v + l)3 3 ^ ^3 3 ^ ^ 



i=l 



32 (2A^3 _^ 2Ar2 - 3A^ - 2) A (-1)^ _ ^ ^}=i 



3Ar(A^ + 1)(A' +2) ^ i2 3^ i + 2 

i=l i=l 



N 

256 



^ . i)-£-.4 ^ 128 (E-.. i-Jii) (E.1. 



3 ^ i 3 ^ i + 2 

1=1 1=1 



Finally, we use J. Ablinger's HarmonicSums package [4-8] , which transforms this expression to 



^The package refers to algorithms and methods from [14,21,29,30,49,50]. 
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the harmonic sum notation: 

64(-l)^ {AN + 1) 2 (270iV^ + SlOiV^ - 46m^ - 1392A^'^ - 211N^ - 20m^ - 156N + 144) 
9{N + 1)^ 277V4(iV + 1)3 

64 , , , /32 , , 16 (lOAr2 + lOAT + 3) \ 32 (lOTV^ + iQAT - 3) 

, / 16 (16A^2 + lOAT - 3) 320 , , 64 , \ 8 (MiV^ + 14iV + 3) , , 
+ W ( 9iV^(iV + ip - + T^^t'" j - 3A^(iV + l) 

1336 ^,,,, 64^ 32(-l)^ \ 32^ , 64^ 128^ 

+ ^52(iV) + y5_2,i(iV) - 3^^^;^3 ) - y52,-2(iV) " y53,i(iV) - — 5_2,l,l(iV). 

We emphasize that the harmonic sums in this expression are algebraically independent. The 

algebraic independence could be accomplished with the Sigma package; out of convenience and 
efficiency we used the HarmonicSums package which contains among various other features the 
harmonic sum relations of [21 J. 

Example 6 The derived sum expression of C^^^^g from Example 3 contains sums of the 
form (22) with depth k — 35. In around four days and 20 hours this expression could be sim- 
plified to an expression in terms of 65 sums that satisfy the properties 1-3 from above. Among 
them there are J^.! sums with depth two; typical examples are 






and 2^^^ ' ■ (24) 

k=l " k=l 

Only one sum of nested depth three has been used, namely 

j 

N 

Z 

k=l 

We emphasize that these sums are constructed in such a way that the difference field algo- 
rithms [20] work most efficiently: The less nested the sums are, the more efficient our algorithms 
work. E.g., if we switch to harm,onic sum notation with Ablinger's HarmonicSum package, the 
first sum in (24) can be rewritten as 

S_siN)S,iN) - S^2,i{N)Ss{N) + S_2{N){lS,{N)^ - IS^iNf - | S,{N)) 

+ 54,-2(iV) - 35_3,L2(iV) - 35_3,2,l(iV) + 35_2 (iV) ^2,1,1 (iV) " ^3,1,-2(A^) 

+ 65_3,i,i,i(^) + 3^-2,i,i,2(iV) + 35_2,i,2,i(iV) - 65_2,i,i,i,i(iV); 

the involved sums have nested depths up to five. With such representations the algorithms in 
Sigma work much slower, or might even fail for our specific input due to time and memory 
limitations. 
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4 3-Loop Anomalous Dimensions and Wilson Coefficients 

In the following we apply the method described in the previous section to unfold all the unpolar- 
ized QCD anomalous dimensions and Wilson coefficients to 3-loop order from a series of Mellin 
moments. This sequence is calculated using the relations given in [3-5] for the different quanti- 
ties per color factor and factors given by ^-values. We will need rather high MeUin moments N. 
The corresponding harmonic sums cannot be calculated by summer [14] directly, but have to be 
evaluated recursively, 

Sa,,a,,...aAN+l) = + I) + Sa.,a,,...a,{N) . (25) 

We used a Maple code for this. The highest moment to be calculated is = 5114 for the Cp- 
contribution to the 3-loop Wilson coefficient C2,q- Its recursive computation requires roughly 
3 GB of memory and 270 min computational time on a 2 GHz processor. It is given by a fraction 
with 13888 numerator and 13881 denominator digits. The set of moments has a size of 69 MB. 
The determination of most of the other inputs sets requires far less resources. 

In Tables 1-3 we summarize the run parameters for the individual color- and (^-contributions 
to the splitting functions and in Tables 4-8 to the Wilson coefficients in unpolarized deeply 
inelastic scattering up to 3-loop order. We specify the number of moments needed on input and 
the order, degree, and length of the recurrence derived. For the solution we compare the number 
of harmonic sums in Refs. [3-5] and in the present calculation. The computation times needed 
to establish and to solve the recurrences are also given. 

To give some example for the rise of complexity for different orders in the coupling constant, 
we compare the contributions to P*-*^^ ^+1. In case of the anomalous dimensions the largest 

amount of moments needed is n = 19 for P„Jr< (n), n = 181 for P^^K,o(n). and n = 1393 for 

99,^ A ^ ' ' 99fiA ^ ' ' 

P^'^^fji (n). The order and degree of the recurrences found are exactly, resp. nearly, the same for 

-^jVS^(^)- ^'^^ non-singlet anomalous dimensions and the singlet anomalous dimensions and 
Pgq^gg{n) ordcr and degree of the difference equation arc larger than in case of Pqg\n). The total 
computation time needed for all anomalous dimensions amount to less than 18 h. The largest 
number of harmonic sums contributing is 26. There are significant reductions in their number 
comparing to the representation given in the attachment to [3, 4] ^ It amounts to a factor of 
two or larger, except in case of the very small recurrences. In the non-singlet case P^g'^s (n) 

the number reduces from 68 to 26. A large reduction is obtained for P^^^'f, in) from 130 to 21 
harmonic sums. 

For the Wilson coefficients C'^2\c'-^^^'^t\c'^ Ca^^^ ^"^^ CfC^ (^) weeks of compu- 
tation time is needed in each case requiring < 10Gb on a 2 GHz processor. The number of 
necessary harmonic sums is 60, reducing from 290 in [5] for C^2\c'^Ca "^^^^ ^'^^ number of 
all harmonic sums not containing the index {—1} up to weight w = 6 after algebraic reduction, 
cf. [30]. 

If one compares the number of harmonic sums obtained in the present calculation after 
the algebraic reduction yields groups characterized by clusters of 58-60, 26-29, 11-15 and cases 
with a number of sums below 10, up to very few exceptions. As this pattern is the same 

^Hcrc, the linear representation given in the text has been reduced already, following an idea of one of the 
present authors. 
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for quite different quantities, it may be related rather to the topology, but the color- or field- 
structure of the respective diagrams. This pattern is not seen counting the harmonic sums in 
the representation of Ref . [3-5] . 

In case of the smaller recurrences the time needed for their derivation is usually shorter than 
that for its solution. Conversely, for the larger recurrences the time required to establish them 
and the solution time behave roughly like 4(3):1. The total computation time amounted to 
110.3 CPU days. Concerning the size of the different problems to be dealt with a naive fivefold 
parallelization was possible. Here wc did not yet consider parallclization w.r.t. the number of 
primes Np chosen, which would significantly reduce the computational time, of the Cp term of 

(3) 

C2q, with Np = 140, discussed above and for other comparably large contributions. 

In course of solving the recurrences we reduce the harmonic sums appearing algebraically, [21] , 
and can express all results in terms of the following harmonics sums: 

^1 

'S'2, S-2 

53, S_3 
S2,l, Si^-2 

54, S-4 

S3,l, S-s^i, S2-2 
S2,1,1,S -2,1,1 

55, S-5 

'S'4,1, 'S'-4,l, 'S'3,-2, 'S'3,2, 'S'-3,2, 'S'-3,-2 

'S'3,1,1, 'S'-3,l,l, 'S'2,2,1, 'S'-2,l,-2, 'S'2,1,-2, 'S'-2,2,l 

'S2,1,1,1,'S'-2,1,1,1 

Se, S-Q 

'S'5,1, 'S'-5,l, 'S'4,2, 'S'4,-2, 'S'-4,2, 'S'-4,-2, 'S'-3,3 

'S4,l,l, 'S_4,l,l, 'S3,2,l, 'S2,3,l, 'S'-3,2,l, 'S'-3,l,2, 'S'-2,3,l, 'S3,l,-2, 'S'-3,l,-2, 'S'-3,-2,l, S-2,2,2, 'S'2,-2,-2 
'S'3,1,1,1, 'S'-3,l,l,l, 'S'2,2,1, 1, 'S'-2,-2,l,l, 'S'2,-2,1,1, 'S_2,2,l,l, 'S'-2,l,l,2 
'S'2,1,1,1,1,'S'-2,1,1,1,1 ■ 

The 3-loop Wilson coefficients require the complete set of possible functions up to w = 6. This 
representation can be further reduced using the structural relations [29, 30] to: 

Si 

S2,l, S-2,1 

S~3,l 

'S'2,1,1, 'S'_2,l,l 
5*4,1, >S'_4,i 

'S3,l,l, 'S'-3,l,l, 'S2,2,l, 'S-2,1, -2, 'S2,l,-2, 'S-2,2,1 
'S2,l,l,l, 'S'-2,l,l,l 

'S-5,1 

'-'4,1,15 '-'-4,1,15 '-'3,2,15 '^2,3,15 '5-3,2,15 '5-3,1,25 '5-2,3,15 '-'3,1,-25 '-'-3,1,-25 '-'-3,-2,15 '-'-2,2,25 '-'2,-2,-2 
'S3 ,1,1,15 'S-3, 1,1,1, .S2,2,l,l5 'S-2,-2,l,l5 'S2,-2,l,l5 'S-2,2,1, 15 'S-2, 1,1,2 
'S2,l, 1,1, 15 'S-2,1, 1,1,1 • 
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In [29,30] we applied a slightly different basis referring to S-2,2,-2 instead of 82-2,-2 and to 5'2,-3,i 
instead of 5'_3^i^2, which is algebraically equivalent. These 38 functions can be represented by 35 
basic Mellin transforms. 

The ab-initio calculation of moments for the quantities considered in the present paper can 
be performed by codes like mincer and MATAD [51] available for physics calculations. Both the 
computational time and memory requests rise drastically going to higher values of N. In case 
of mincer both parameters increase by a factor of ~ 5 enlarging N ^ N + 2. Comparable, but 
slightly larger factors are obtained for MATAD. In the well-known leading order case, enough 
moments may be provided for our procedure. Already for some color projections of the next- 
to-leading order corrections, this is no longer the case, [53], since around 150 initial values are 
needed. For the 3-loop anomalous dimensions and Wilson coefficients N — 16 can be reached 
with computation times of the order of 0.5-1 CPU year, cf. [9]. The codes [51] still may be 
improved. However, the power-growth going to higher moments will basically remain due to the 
algorithms used. The method presented in this paper can therefore not be applied to wliole color- 
factor contributions for the anomalous dimensions and Wilson coefficients at the 3-loop level. 
They may, however, be useful in solving medium-size problems. In view of constructing general 
methods suitable to evaluate single scale quantities, methods to evaluate the fixed moments for 
these quantities at far lower expenses have to be developed. 

To illustrate the results of the present calculation, the non-singlet anomalous dimensions to 
O(a^) are given as an example in the appendix. The relations for all unpolarized anomalous 
dimensions and Wilson coefficients, separated according to the corresponding color- and C- value 
terms, are attached to this paper in FORM- and Mathematica ffies. The FORM-codes provide a 
check of our relations with the moments calculated in Ref . [7, 8] . 

5 Conclusions 

We established a general algorithm to calculate the exact expression for single scale quantities 
from a finite, suitably large number of moments, which are zero scale quantities. The latter 
ones are much more easily calculable than single scale quantities. We applied the method to 
the anomalous dimensions and Wilson coefficients up to 3-loop order. Hereby we compactified 
their representation exploiting all algebraic relations between the harmonic sums. The 3-loop 
Wilson coefficients require the whole set of basic harmonic sums in the sub-algebra spanned by 
the index set to w = 6 without i = —1. A further compactification can be obtained using the 
structural relations between the harmonic sums. After algebraic reduction the number of the 
harmonic sums contributing clusters in several classes mainly determined by the topology of the 
graphs and widely independent of the color- and field structure of the respective contributions. 
The CPU time for the whole problem amounted to about four months using 2 GHz processors 
and < 10 GB of memory were needed. The problem can be naively parallelized fivefold. The 
real computational time needed to establish the recurrences can be shortened further running 
Chinese remaindering in parallel. 

To solve 3-loop problems for whole color factor contributions is not possible at present, since 
the number of required moments is too large for the methods available. Methods to evaluate the 
fixed moments for these quantities to high order at far lower expenses have still to be developed. 

We established and solved the recurrences for all color resp. C-pi'oj actions at once, which 
forms a rather voluminous problem. Yet we showed that rather large difference equations [order 
35; degree ~ 1000], which occur for the most advanced problems in Quantum Field Theory, can 
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be reliably and fast established and solved unconditionally. 
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Table 1: Run parameters for the unfolding of the non-singlet anomalous dimensions 





number of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a [b] 




Pns,o 


14 


2 


3 


0.05 


0.087 


1 [1] 


0.55 


P~ 

NS.l.C]^ 


142 


5 


31 


3.32 


4.666 


6 [10] 


7.45 


^NS.I.C'aC'f 


109 


4 


24 


1.91 


2.834 


6 [7] 


6.28 


A .SM.C'fAV 


24 


2 




O.LS 


0.271 


2 [2] 


0.92 


P+ 
JV5,1,C|, 


142 


5 


31 


3.35 


4.707 


6 [10] 


7.45 


P+ 


109 


4 


23 


1.88 


2.703 


6 [7] 


6.23 


P+ 

^NS,1,CfNf 


24 


2 


7 


0.09 


0.271 


2 [2] 


0.89 


P~ 
NS,2,C% 


1079 


16 


192 


3152.19 


529.802 


26 [68] 


1194.41 




48 


3 


11 


0.49 


0.643 


1 [1] 


1.56 


P~ 
NS,2,CaCI. 


974 


15 


181 


1736.08 


450.919 


26 [62] 


1194.41 


P~ 

^ NS,2,CaCIC3 


48 


3 


11 


0.53 


0.643 


1 [1] 


1.53 


P~ 

NS.2,C^Cf 


749 


12 


147 


1004.12 


242.892 


26 [62] 


1100.88 


NS.2,C\CFi3 


48 


3 


11 


0.56 


0.643 




1.56 


NS.2,CfNI 


39 


2 


11 


0.31 


0.369 


3 [3] 


1.20 


P~ 

NS,2,C%Nf 


377 


8 


68 


76.34 


33.946 


11 [24] 


72.22 


^NS,2,Cj,NFC3 


14 


2 


3 


0.12 


0.101 




0.53 


PnS,2,CaCfNf 


356 


7 


62 


65.25 


23.830 


11 [20] 


52.67 


PnS,2,CaCfNf(3 


14 


2 


3 


0.12 


0.101 




0.55 


P+ 

^NS,2,C% 


1079 


16 


192 


4713.27 


527.094 


26 [68] 


1165.22 


P+ 

^NS,2,CI,(3 


48 


3 


11 


0.55 


0.643 




1.562 


P+ 
NS,2,CaC^ 


974 


15 


178 


1715.03 


442.031 


26[62] 


889.047 


P+ 

^NS,2,CaCI,C3 


48 


3 


11 


0.61 


0.643 




1.531 


P+ 


749 


12 


146 


991.22 


240.325 


26[50] 


516.812 


P+ 

^NS,2,CICfC3 


48 


3 


11 


0.61 


0.643 




1.593 


P+ 


377 


8 


69 


111.38 


33.872 


11 [24] 


71.235 


P+ 


14 


2 


3 


0.15 


0.101 




0.531 


P+ 


307 


7 


61 


48.62 


23.808 


11 [24] 


71.235 


P+ 


14 


2 


3 


0.15 


0.101 




0.547 


P+ 

Ar5,2,CFAf^ 


39 


2 


11 


0.40 


0.369 


3 [3] 


1.172 




459 


7 


87 


239.62 


0.369 


5 [20] 


32.5 
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Table 2: Run parameters for the unfolding of the singlet anomalous dimensions 





number of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a[b] 




r>PS 

1 ATiT-r^c 

X , J V P ±r 


24 


1 


8 


0.19 


0.204 


0[0] 


0.244 


pPS 

Z,i V O F 


109 


3 


26 


6.32 


1.891 


2 [8] 


2.812 


pPS 

-t O AT /~< 

2 , iV F O A F 


566 


9 


115 


425.44 


100.414 


7 [40] 

L J 


111.52 




19 


1 


6 


0.42 


0.117 


0[0] 


0.204 


pPS 


237 


5 


50 


32.75 


0.117 


4[24] 


14.601 


pPS 


19 


1 


6 


0.41 


12.163 


0[0] 


0.200 


P n 


11 


1 


3 


0.02 


0.061 


[0] 


0.16 






A 


9Q 


O.Ol 


0. ( uy 


[81 
[J 


4 879 
4.0 i z 


Pqg,l,NpCF 


uo 


o 




U.wO 


n Q^l 


L J 


9 nns 


P 

^qg,2,NlCA 




D 


04 


40. oU 


1 7 A(\'\ 
1 i .4Uo 


D[Z4J 


91 00*? 

zi .yyo 


^qg.2,NlCF 


1 0/1 





A 1 
41 


97 7Pi 
z ( . ( 


7 01 1 

/ .yii 


Q [1 ^1 
o[10J 


8 n9i 

o.Uzl 


P 






9ni 

ZUl 


00Z1.40 


cry cqc 

00 / .000 


1 *? f88l 
10 [ooj 


8/18 8^ 


P o 




9 


1 1 


u.ou 


n ADS 
U.4UO 


1 ["^1 
1[0J 


n Q'?9 

u.yoz 


Pqg,2,NFCACF 




1 

±o 




'^Qfi'? fi9 


'i'i9 inn 

00^. XULJ 


1 9 rs4i 

x^ l^^l 


71 4 4^ 


Pqg,2,NFCACFC3 




9 


1 1 

X X 


1 09 
± .uz 


n 400 


1 \6\ 


u.yo 


Pqg,2,NFCl 


849 


12 


143 


1 337 36 


261 804 


13 [66] 


387.6 


Pqg,2,NFC],C3 


17 


1 


5 


0.29 


0.093 


[0] 


0.15 


P „ 


1 1 

1 1 


1 
1 


•3 
O 


U.UO 


U.UOZ 


U [UJ 


n 1 ^ 
u. 10 






o 




fl 79 

U. ( z 




L J 


1 Q94 


p , ^ ^ 


IZO 


/I 


Ol 


4.00 


4 nl^o 
4.uoy 


■^ri 9i 

0[1ZJ 


O.UOO 


J </(/.l,Aj.-('j.- 




2 






V ' ■ X . / z 






P 

^gq,2,Cj, 


/Uo 


1 1 
i i 


1 14t 


097 89 
yz / .oZ 


1 K9 '?9n 
IDZ.oZU 


io[oyj 


Z40.0O 






9 


Q 


n 7Q 


n 981 

u. ^ox 


i[dj 


77fi 


P 


iUoo 


10 


ZUr3 


0097 09 
OOZ ( .oZ 


0OO.O4D 


1 9 fO'^l 

iz [yoj 


8'^n 9 
ooU.Z 


P 




9 


Q 

y 


n 81 

U.Ol 


n 981 
U.ZoX 


1 ['^1 
i[dj 


n 77R 


P 

-f^g,2,C^C|, 


1 nc7 
iUo ( 


10 


1 OQ 

lyo 


olo4.1o 


COQ C07 

Ozo.oz I 


1 A r7f;1 
14 [(Oj 


000. ol 


P 

^gq.2,CAClCi 


OO 


9 


Q 

y 


U.OD 


n 981 


1 ["^1 


n 77R 


Pgq,2,NFCj, 


OOt? 


7 


fiQ 




•^n fi9fi 


^^[9111 
1^01 


oo.oou 


Pgq,2,NFClC,3 


1 1 

± ± 


1 
± 


Q 
O 


n 94 


n nfi9 


nrni 

L J 


n 1 f^9 




1087 


15 


194 


3201.52 


58.943 


17[87] 


714.51 


Pgq,2,NFCACFC3 


11 


1 


3 


0.21 


0.062 


0[0] 


0.156 


Pgq,2,NlCF 


41 


3 


9 


1.04 


0.445 


2 [6] 


1.216 


Pgg,o 


19 


2 


5 


0.04 


0.166 


1 [1] 


0.65 


Pgg,i,cl 


181 


5 


45 


12.07 


9.053 


6 [17] 


11.62 


Pgg,l,NFCA 


29 


2 


9 


0.23 


0.395 




0.856 


Pgg,l,NFCF 


31 


1 


11 


0.23 


0.228 


0[0] 


0.240 
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Table 3: Run parameters for the unfolding of the singlet anomalous dimensions (continued) 





number of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a [b] 






1393 


16 


277 


12432.80 


1087.615 


21 [130] 


2419.04 


P„„ 2 Ar„(72 


439 


8 


88 


237.82 


57.291 


7 [35] 


55.52 


PggXNpClCs 


15 


1 


4 


0.31 


0.073 


0[0] 


0.156 


Pgg,2,NFC\ 


782 


11 


148 


1638.62 


205.980 


6 [31] 


160.89 


Pgg,2,NFC\ii3 


29 


2 


9 


0.66 


0.308 


1[1] 


0.796 


Pgg,2NFCACF 


749 


10 


127 


1169.37 


146.921 


7 [40] 


128.37 


PggaNFCACFCi 


29 


2 


9 


0.72 


0.305 




0.828 


Pgg,2,N%CA 


55 


2 


17 


4.53 


0.979 


1[4] 


1.092 


Pgg,2,N%CF 


109 


3 


26 


6.74 


2.483 


2[12] 


2.668 



Table 4: Run parameters for the unfolding of the unpolarized pure-singlet Wilson Coefficients 





number of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a[b] 




^2,PS,CfNf 


209 


5 


42 


20.85 


8.422 


3[16] 


7.70 


^(3) 

^2,PS,C%Nf 
^(3) 

'~'2,PS,Cj.NFC3 
^(3) 

2,PS,CI,NfU 

^2,PS,CfCaNf 

^2,PS,CfCaNfC3 
^(3) 

^2,PS,CfCaNf<:4. 

^2,PS.CfNI 
^(3) 

^2,PS,CfNIC3 


1847 
65 
19 

2023 
71 
19 
479 
19 


19 
2 
1 

20 
2 
1 
8 
1 


334 
20 
6 

368 
21 
6 

103 
6 


41001.00 
0.69 
0.08 

54873.80 
0.82 
0.08 
629.05 
0.09 


1989.043 
1.124 
0.117 

2670.459 
1.429 
0.117 
75.646 
0.122 


14 [122] 
1 [6] 
0[0] 

14 [126] 
1[6] 
0[0] 
5 [34] 
0[0] 


2701.04 
0.92 
0.14 

4589.31 
0.97 
0.14 
53.28 
0.14 


^(2) 

^L.PS.CfNf 


41 


2 


11 


0.20 


0.384 


1[4] 


0.88 


^L,PS,C1.Nf 
^L,PS,C],NfC3 

'~'L,PS,CfCaNf 
ri(3) 

'~'L,PS,CfCaNfC3 
^(3) 

^L,PS,CfN^ 


869 
35 

840 
35 

224 


11 
2 

11 
2 
5 


162 
10 

153 
10 
52 


4411.20 
0.17 

2005.44 
0.17 
72.64 


250.352 
0.406 

231.837 
0.403 
12.440 


8 [62] 
1 [5] 

8 [64] 
1[5] 

3 [13] 


163.17 
0.63 

153.99 
0.59 
9.87 
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Table 5: Run parameters for the unfolding of the unpolarized quarkonic Wilson Coefficients for the 
structure function F2{x,Q'^). 





nuinbcr of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a[b] 




^(1) 


31 


3 


6 


0.26 


0.429 


2 [3] 


1.47 




689 


11 


137 


1134.10 


177.806 


13[39] 


258.24 


^2,q,CACF 

C2,q,CACF(3 

C2,q,NFCF 


15 
545 
15 
71 


2 
10 
2 
4 


3 
121 
3 

16 


0.27 
413.33 
0.27 

2.68 


0.100 
127.893 
0.112 

1.655 


12[35] 

4[10] 


0.54 
178.73 
0.55 

3.95 


^(3) 


5114 


35 


938 


1.79 xlO^ 


30394.173 


58[289] 


509242 


r-(3) 

^2,g,C|,f4 
ri(3) 

^2,g,C|,C5 
/^(3) 

'-'2,qr,C|.CA 
^(3) 

'-^2,q,C|.CAC3 
r^(3) 

'-^2,q,C|.CAC4 
^(3) 

'-'2,9,C|,CaC5 
/^(3) 


284 
65 
19 
5059 
284 
48 
19 
4564 


8 
3 
2 

35 
8 
3 
2 

33 


64 
11 

5 

930 
64 
11 
5 

863 


31.02 
2.62 
0.08 
1.69 xlO^ 
34.00 
0.32 
0.08 
1.39 xlO'^ 


32.363 
0.163 
0.163 
30122.380 
33.400 
0.643 
0.167 
24567.518 


6 [18] 

1 [1] 
1 [1] 
60 [290] 

7 [18] 
60 [258] 


27.60 
1.47 
0.47 
0.478 xlO^ 
28.53 
1.01 
0.42 
0.349 xlO^ 


/^(3) 

'^2,7,CfC2C3 
/^(3) 

rm 

^2,q,CFClQ, 
/^(3) 

^2,7,C|,JVf 


284 
48 
19 
1762 


8 
3 
2 
20 


63 
11 
5 

348 


26.83 
0.32 
0.08 
40237.45 


29.918 
0.643 
0.175 
2339.516 


7 [17] 
28 [107] 


30.46 
1.01 
0.42 
7548.56 


/^(3) 

2,q,CpN fC,3 

^2,q,ClNFU 
/^(3) 

^2,(},CFCAAfF 

/^(3) 

^2,g,CFCAAfFC3 

ri(3) 

'-^2,g,CFCAAfFC4 
^(3) 

^2,q.CFNl 
^2,q,CFNlC:i 


87 
15 
1847 
89 
15 
131 
15 


4 
2 
20 
4 
2 
5 
2 


21 
3 

360 
24 
3 
30 

3 


1.94 
0.07 
47661.64 

2.47 
0.06 
58.00 
0.06 


2.354 
0.101 
2507.362 

2.935 
0.101 
5.347 
0.101 


3 [5] 

28 [111] 

3 [8] 

7 [22] 
1 [1] 


2.83 
0.34 
7525.89 

3.19 
0.34 
12.22 
0.38 


^(3) 

^2,q,dabcC,3 
^(3) 

"-^2,g,do6cC5 


1199 
109 
8 


14 
4 
1 


242 
25 
2 


6583.27 
2.33 
0.03 


738.498 
3.164 
0.041 


15 [62] 
2 [7] 
0[0] 


841.24 
2.40 
0.10 
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Table 6: Run parameters for the unfolding of the unpolarized quarkonic Wilson Coefficients for the 
structure function Fl{x,Q'^). 





number of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a[b] 




-o(l) 

L q Cf 


5 


1 


1 


0.02 


0.033 


0[0] 


0.12 


L,q,Cl 


203 


5 


51 


7.86 


12.381 


6[11] 


17.21 


L,q,ClC3 


5 


1 


1 


0.02 


0.033 


0[0] 


0/13 


^L.q.CACp 


159 


4 


43 


4.32 


7.624 


5[8] 


11.43 


^L,q,CACF(i3 


5 


1 


1 


0.02 


0.033 


0[0] 


0.12 


'-'L,q,CpNp 


19 


2 


4 


0.05 


0.134 


1[6] 


5.30 




2419 


22 


472 


110504.41 


4555.679 


27 [142] 


19060.10 




131 


5 


34 


3.52 


6.257 


3 [10] 


7.40 


^(3) 


11 


1 


3 


0.05 


0.069 


0[0] 


0.14 


p(3) 
L a CiCA 

'/ y '/ ^ ^ 


2551 


23 


486 


124064.39 


5176.054 


27 [144] 


24614.00 




131 


5 


34 


4.51 


6.807 


3 [10] 


7.39 


^(3) 


11 


1 


3 


0.05 


0.069 


0[0] 


0.14 


^(3) 

'-"L^q^CpC^ 


1803 


18 


344 


42500.82 


2064.227 


27 [109] 


6269.33 


/^(3) 

^L,q,CFClC3 


131 


5 


31 


3.50 


5.463 


2 [10] 


6.32 


-o(3) 

'^L,q,CpCl(^ 


11 


1 


3 


0.05 


0.069 


0[0] 


0.15 


^L,q,CFCANF 


1014 


14 


203 


4041.82 


539.901 


13 [58] 


896.70 


ri(3) 

'-'L,q,CFCANF(i3 


41 


2 


12 


0.19 


0.518 


1[5] 


0.92 


^(3) 

^L,q,Cl.NF 


959 


13 


188 


3507.92 


400.784 


13 [51] 


769.90 


^(3) 

'-'L,q,C%NFC3 


29 


2 


8 


0.15 


6.257 




0.85 


(-1(3) 


47 


3 


10 


1.58 


0.498 


2 [4] 


1.45 


r'(3) 

^L,q,dabcNF 


989 


12 


184 


3536.04 


371.269 


15[60] 


384.00 


MS) 

L .q.dahcN F 


89 


4 


18 


1.90 


2.034 


2 [7] 


2.68 


p(3)^ 


5 


1 


1 


0.02 


0.033 


0[0] 


0.12 
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Table 7: Run parameters for the unfolding of the unpolarized gluonic Wilson Coefficients for the 
structure function F2{x,Q'^). 





number of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a[b] 






24 


2 


6 


0.14 


0.191 


1[3] 


0.72 


^(2) 
^(2) 

'-'2,g,NpCAC3 
^2.g.NFCF 


459 
8 

419 

8 


9 
1 
8 
1 


93 
2 

91 
2 


202.96 
0.11 

207.98 
0.14 


73.022 
0.038 

63.468 
0.038 


7[35] 
0[0] 

7 [32] 
0[0] 


70.77 
0.14 

59.78 
0.14 


^2,g,C%NF 


3464 


28 


658 


542132.00 


11742.788 


29 [228] 


65721.40 


^(3) 

M3) 

'''2,g,ClNF<5 
^(3) 

'-'2,g,ClNF 


181 
17 
11 
4014 


6 
1 
1 
30 


42 
5 
3 

739 


25.30 
0.23 
0.20 
869580.00 


12.171 
0.093 
0.067 
16320.095 


3 [14] 
[0] 
[0] 
28 [261] 


7.67 
0.12 
0.15 
97289.30 


r(3) 

^2,g,C\NpC3 
^2,g,C\NFC4. 

r<(3) 
r(3) 

'^2,g,CFCANF 
^(3) 

'^2,g,CFCANFC3 

r(3) 

^2,g,CFCANFCi 
-o(3) 

^2,g,CFCANFia 

r(3) 

'^2,g,CFNl 


194 
39 
11 

4014 
194 
39 
11 

1553 


6 
2 
1 

30 
6 
2 
1 

16 


44 
11 

3 

747 
43 
11 
3 

285 


42.39 
0.74 
0.17 
889246.00 
41.81 
0.61 
0.17 
22235.00 


13.263 
0.408 
0.063 
16640.997 
12.999 
0.409 
0.068 
1181.805 


3 [15] 
1[3] 
[0] 
29 [264] 
3 [15] 
1[3] 
[0] 
13 [101] 


8.01 
0.67 
0.13 
100830.00 
7.90 
0.66 
0.11 
1506.79 


^-(3) 

r^(3) 

'-'2,g,CAN^ 
^(3) 

'-'2,g,CAN^C3 


55 
1329 
39 


2 

17 
2 


16 
259 
11 


2.81 
10692.80 
2.48 


0.962 
1033.138 
0.666 


1[3] 
13 [96] 
1[3] 


0.70 
1162.99 
0.70 


r^(3) 

^2,g,dabcNF 

r<(3) 

^2.g.dabcNfC3 
^(3) 

^2,g,dabcNFC5 


1403 
142 
19 


15 
5 
1 


282 
37 
7 


13951.90 

8.54 
0.30 


1048.336 
7.177 
0.139 


19 [81] 

2 [12] 
[0] 


2668.66 
6.74 
0.14 
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Tabic 8: Run parameters for the unfolding of the unpolarized gluonic Wilson Coefficients for the 
structure function Fl{x,Q'^). 





number of 


order of 


degree of 


total time 


length of 


number of 


solution 




terms 


recurrence 


recurrence 


[sec] 


recurrence 


harm, sums 


time [sec] 




needed 








[kbyte] 


a[b] 






5 


1 


1 


0.02 


0.033 


0[0] 


0.13 


L,g. C p Mp 
^L,g,CANF 


153 

109 


4 

4 


38 

25 


4.15 
1.31 


5.941 
2.731 


2 [6] 
3[10] 


5.30 
4.22 


^(3) 

^L,g,ClNp 


1679 


17 


314 


48496.50 


1498.918 


16 [100] 


2019.46 


^(3) 
^(3) 

^L,g,ClNi. 


120 

5 

1671 


4 
1 
17 


28 
1 

302 


3.64 
0.02 
29219.30 


3.967 
0.033 
1392.205 


2 [8] 
0[0] 
16 [112] 


3.23 
0.09 
2012.38 


'-'L,g,ClNpC3 

^L,g,ClNFQ. 

^L,g,CFCANF 
^(3) 

'^L,g,CFCANpC3 
^(3) 

^L,g,CFCANF(:5 
^(3) 

'^L,g,CFNl 

^L,g,CFNlC.3 

^L,g,GANl 
^(3) 


109 
5 

1935 
120 
5 

699 
15 

419 

5 


4 
1 
18 
4 
1 
9 
1 
8 
1 


24 
1 

351 
28 
1 

140 
4 
90 
1 


2.46 
0.03 
44671.90 
4.43 
0.03 
1350.09 
0.06 
526.25 
0.02 


3.007 
0.033 
2036.550 
4.154 
0.033 
140.949 
0.074 
57.569 
0.033 


2[8] 
0[0] 
16 [116] 
2 [8] 
0[0] 
6 [35] 
0[0] 
6 [30] 
0[0] 


2.836 
0.11 
3510.31 
3.10 
0.11 
108.69 
0.17 
47.40 
0.08 


^(3) 

'~^L,g,dabcNF 

L,g,dabcNF (^3 

n{3) 

L,g,dabcNFC5 


1109 
129 
11 


13 
5 
2 


231 
27 
2 


10155.40 
2.18 
0.06 


618.402 
3.858 
0.074 


18 [75] 
2 [11] 
0[0] 


1714.70 
4.09 
0.12 
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6 Appendix: The non-singlet anomalous dimensions 

The non-singlet anomalous dimensions Pgq^ and Pq'^(n)\k=i,2 are given by 



Pi fin) 



P^fin) 



4^1 



?>r? + 3n + 2 



+ 1) 

3n'"' + 9//''' + 9// ^ - 5//-^ - 2 1/r - 32// - 24 



(26) 



2n3(n + l)3 



- 165_3 



1(3 



n(n-l-l) / yn^yn^iy J n{n-\-l) 

16(-1)""' 



-16^3 + 325_2,i + 
+ CaCf 



hlrv' + 102n^ + 65571^ + 48471^ + 12n + 144 



18n3(n + l)2 



+-5_2 ( 16>5i 



44 8(-l)" 

( , ^\ I q"'^2 + 8S'3 - 165'_2,i - 7 — TT 

n(n + 1)/ 3 (n + 1)- 



3n^ + 6n3 + A7n'^ + 20n - 12 40 8 

9n2(n + l)2 y ' + 3 ' 



+5- 



+ 9n^ + + 59n^ + 40/1^ + 32n + 8 



2n3(n + l)3 

n(n-l-l) / yn^'^n + 1)^ 

16(-1)' 



16S. 
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+ ^1 



3n(n + 1) 

2 (245n^ + 980n^ + 1542/1^ + 1524n^ + 851n^ + 100^3 + 3671^ + 22n - 6) 

3n^{n + 1)4 



8(lln^ + lln-8),g3 _ _ _ 32 (lln^ + lln - 24) ,g_2,i 

n{n + l) ^ 3n(n + l) 

+ 51252,-2 + 256^3,1 + 10245_2,i,i ] + 7685_2,2,i + 768^2,1,-2, - 384^3,1,1 

/ 28 



1536S'_2,i,i,i 



+ 



+ 



+ 



12 (5n^ + 10^3 + - An - A) 
v?{n + 1)^ 



+ CfNI 



+ ClNp 



Sln^ + 153n^ + 57n^ + 35n3 + 96n2 + 16n - 24 16 80 16 

27n3(n+l)3 27^^ " 27^^ + 9 ^ 

_ 32 2 _ 4 (15n^ + 30^3 + 7971^ + 16n - 24) ^2 
2 9n2(n + l)2 



207n^ + 828n^ + 144371^ + 112371^ - 38n^ - 77971^ - 63271^ + 120 128 



+ 5-3 
+ S^2 
+ 5l 



97z4(rz + ly 
32(10„' + 10n + 3)_g4 X 

9n(n + l) 3 7 

32 (16n2 + lOn - 3) 640 _ 12 



64^1 128(4rz + l) 



3(rz+l)3 9(n+l)^ 



9n2(n + 1)2 



+ 



9 



-Si 



^2 + 



16 (29ri2 + 29n + 12) ^3 _ 128 , 
9n(nTT) 3~' 



2 (165/1^ + 33072^ + 165n3 + 160rz2 - 1671 - 96) 320 128 128 
64 (107z2 + 1071 - 3) 5_2 1 64 64 256 



32^1 - 



8 (37z2 + 372 + 2) 



+ CaCfNf 



64 



n{n + 1) 

2 (270 + 81071^ - 46371^ - 139271^ - 21171^ - 20671^ - 15671 + 144) 



277i4(n + l)3 



32 16 (1071^ + 1071 + 3) \ 



1336 



S2 + 'S'_2 



3"' 972(7Z+1) 

16 (167i2 + 107Z - 3) 320 



) 



+ (-!)■' 



/64(47z + l) 32^1 



V 9(72 +1)4 3(72+1)3 



32 (1072^ + 107Z - 3) S-2,1 



9n{n + 1) 

64 

+ 1653+y5_2,l 



64 \ _ 8(1472^ + 1472 + 3)^3 80 
9^37 372(72 + 1) 3 ^ 

4 (20972^ + 62772^ + 62772^ + 281723 + 3672^ + 3672 + 18) 



27723(72 + 1) = 



32 64 128 

- -2 - y^3,i — ^<^_2,i,i 



(30) 



C'f 



+ s. 



64 



72(72 + 1) 
64 (372^ + 772 + 5) 

722(72+1)2 



- 128^1 1 S^_2 + 



16 (3//^^ + 9/;'' + + ir' + 2ir + 1// + 2) 



1408^2 



723(72 + 1)3 

64 (372^ + 372 - 11) ^2 
72(72+1) 



+ 1536^3 + 1285_2 1 
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' n{n + l) 2n^(n + l)^ 



n(n + 1) / \ n(n + 1) 

+ ^''<^"; + ^" + ^) - 320S_. - 2240S.'| + (-1)" ( "'(^"^ + ""+^' + ^^SS-, 



96(5n + 3)^1 _ 64^2 A 4 (13n^ + 2671^ + 13n^ - 16n - 20) ^3 
^ (n + l)4 + J ^ 

n[n + Ij n(n + 1) 

2/ 32 (3n^ + 3n + l) \ 32 (3n^ - n + 4) 5-2,i 

+ '^^^ n3(n + l)3 n2(n + l)2 

, _ /2(3n6 + 9n5 + 9n^ + 83n3 + 76n2 + 60n+16) , , 

+ -^2 T} —77^ r D4^3 + ^i7D<b_2,l 
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8 (37z2 + 3n + 2) 
n{n + 1) 



325i 



32 
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3n5(n + l)5(n + 2)^ 



+ 



4 (n^ + n + 2) S. 



P9{n)Si 



n2 n + 1 



+ S. 



+ (-1) 



8^1 (n^ + n + 2)^ 



4 (n^ + 3n^ 



3n4(n + l)4(n + 2)3 
21n'^ - 23n2 - 12n - 4) 



n 
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16 {5n^ + 29n^ + 78n* + llSn^ + lUn^ + 72n + 16) 5i 



3(n- l)n2(n+ l)3(n + 2)3 

4 (13n^ + 74n^ + ITQn^ + 314/1^ + 644n^ + lOOOn^ + 81671^ + 352n + 64) 
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+ 1)^ 
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(31) 



For brevity we abbreviated Sg{n) = Sg- Here, Ca = N^, Cp = {N^ - l)/(2iVc) are SU{Nc) color 
factors, Np denotes the number of quark flavors and Nc is the number of colors, with iVc = 3 for 
Quantum Chromodynamics. We have accounted for the color factor Tr — 1/2 explicitly, which 
is the same for all groups SU{Nc). dabc denotes a SU{N^ structure constant and the Einstein 
convention is applied calculating dahcd""^^- 

The functions Pi{n) which appear in Eqs. (26-31) are given by 



Pi(n 
A(n 

P^n 

P^{n 

Pe{n 

P,{n 
Ps{n 

PM 



= 29n^° + U^n^ + 130n^ - 146n^ - 479n^ - lln^ - 464n^ - 1748ri=^ - 1600n^ 

-752n - 16 (32) 
= 1359n^° + 679571^^ + 15246n^ + 15646n^ + 3851n^ - 35089n^ - 34648n^ 

+12280n^ + 32592n^ + 17616n + 3456 (33) 
= 4971n^" + 24855n'' + 11770n^ - 86322n^ - 150929n^ - 135893n^ - 85692n^ 
+ -18992n^ + 2282471^ + 1584071 + 259 (34) 
= 2971^° + 14571*^ + 22671^ + 11071^ + 35371^ + 50l7i^ + 97671^ + 9407i^ + 57671^ 

+20871 + 32 (35) 
= 135971^ + 543671^ + 827471^ + 2452471^ + 1110371^ + 1252871^ + 4 12071^ 

-256071 - 1584 (36) 
= 497171^ + 2485571^ + 1076271^ - 5713871^ - 9203371^ - 4090l7i^ + 1069271^ 

+121671 - 2904 (37) 
= 20771^ + 82871^ + 149171^ + 229l7i^ + 133871^ + 45371^ - 87i^ - 16071 - 72 (38) 
= 4(1371^° + 9771*^ + 32671^ + 72071^ + 139971^ + 241 671^ + 301771^ + 241271^ 

+1184rz2 + 336n + 48) (39) 
= 2(971^ + 4l7i^ - 371^ - 50571^ - 171971^ - 295l7i^ - 309271^ - 203271^ 

-76871 - 144) (40) 
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